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Abstract 

We show that the recent measurements of Bs — Bg mass difference, Am^, by D0 and 
CDF cohaborations give very strong constraints on MSSM scenario with large flavor 
mixing in the LL and/or RR sector of down- type squark mass squared matrix. In 
particular, the region with large mixing angle and large mass difference between scalar 
strange and scalar bottom is ruled out by giving too large Anig. The allowed region 
is sensitive to the CP violating phases Si^i^jiy The Amg constraint is most stringent 
on the scenario with both LL and RR mixing. We also predict the time-dependent 
CP asymmetry in Bg ip4> decay and semileptonic asymmetry in Bg iX decay. 
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1 Introduction 



The flavor changing processes in the s — b sector are sensitive probe of new physics (NP) 
beyond the standard model (SM) because they are experimentally the least constrained. 
In the minimal supersymmetric standard model (MSSM), however, the flavor mixing in the 
chirality flipping down-type squarks, s^l{r) — bR(^L), is already strongly constrained by the 
measurement of BR(B Xg'-f). On the other hand, large flavor mixing in the chirality 
conserving ^l{r) — ^l{r) has been largely allowed. Especially the large mixing scenario in 
the 'sr — hn sector has been drawing much interest because it is well motivated by the 
measurement large neutrino mixing and the idea of grand unification [1]. 

Recently D0 and CDF collaborations at Fermilab Tevatron reported the results on the 

measurements of Bg — Bg mass difference [2, 3] 

17 ps"^ < Am, < 21 ps-^ (90% CL), 

Am, = n.SStoil ± 0.07 ps-\ (1) 

respectively. These measured values are consistent with the SM predictions [4, 5] 

Am^^(UTfit) = 21.5 ± 2.6 ps-\ Am^^(CKMfit) = 21.7^1^ ps"^ (2) 

which are obtained from global fits, although the experimental measurements in (1) are 
slightly lower. The implications of Am, measurements have already been considered in 
model independent approach [6, 7, 8], MSSM models [9, 10], Z'-models [11], etc. 

In this paper, we consider the implications of (1) on an MSSM scenario with large 
mixing in the LL and/or RR sector. We do not consider flavor mixing in the LR(RL) 
sector because they are i) are already strongly constrained by BR{B — > Xg^) [12] and 
ii) therefore relatively insensitive to Bg — Bg mixing. We neglect mixing between the 1st 
and 2nd generations which are tightly constrained by K meson decays and K — K mixing, 
and mixing between the 1st and 3rd generations which is also known to be small by the 
measurement of B^ — B^ mixing. 

The paper is organized as follows. In Section 2, the relevant formulas for Bg — Bg mixing 
are presented. In Section 3 we perform numerical analysis and show the constraints imposed 
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on our scenario. With these constraints, in Section 4, we predict the time-dependent CP 
asymmetry in Bs ip(f) decay and the semileptonic asymmetry in Bg —>■ iX decay. We 
conclude in Section 5. 



2 Bs — Bs mixing in the MSSM scenario with large 
LL/RR mixing 



According to the description of our model in Section 1, the scalar down- type mass squared 
matrix in the basis where down quark mass matrix is diagonal is given by [13, 14] 
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/ ~d,2 

' m 



M~ 

d,LL 









-d,2 





'd,2 



1^22 ^23 



V 



0-^(^,2 --^ d.l 



d,LR{RL) 



The Mf can be obtained from ^ by exchanging L 



^ R. We note that this kind 
of scenario is orthogonal to the one with flavor violation controlled only by CKM matrix 
(minimal flavor violation model [15, 8] or the effective SUSY model considered in [16]), 
where large flavor violation in s — 6 is impossible a priori. 



The mass matrix M~ ^ ^ can be diagonalized by 



(4) 



with 









cos^^L sin^^i e 
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(5) 



y — sin^^e '"^ cosUl j 

Similarly, the exchange L -^-^ i? in (5) gives Fr. We restrict —45° < Ql{r) < 45° so that the 
mass eigenstate 's{b) has more strange (beauty) flavor than beauty (strange) flavor. 
The most general effective Hamiltonian for Bg — Bg mixing 
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-f^eflf = ^ CiOi + ^ CiOi 



=1 
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has 8 independent operators 



Ui = 












04 = 








;=i,"3 — 


Oi=l,...3 \l*^R ■ 



Oi=l,...3 — Oj=i^...3 \l*^r . (7) 



The Wilson coefficients for these AS = I^S — 2 operators can be obtained by calculat- 
ing the gluino mediated box diagrams. Since the chargino and neutralino exchanged box 
diagrams are suppressed by the small gauge coupling constants, we neglect them. In the 
scenario we are considering, when we consider only LL (RR) mixing, the SUSY box diagram 
contributes only to Ci (Ci). When both LL and RR mixing exist simultaneously, there are 
also contributions to C4 and C5. However, 2 or (7^3 are not generated at all. Note that 
the induced LR (RL) mixing [17] does not occur, cither, because wc set M~ = Osxs- 

Otherwise, the SUSY parameter space is further constrained depending on tan/3 [17]. The 
analytic formulas for the Wilson coefficients at the MSSM scale are given by 

^M^r - ^ sin 2^^e^(^-+^-) (^(5) (x^,,,, x~,J - ^(5) {x~,^~^, x,,,,) 

~f4:(5){XsR,g,Xl^~) + f4(5){Xsji,g,Xsj^^g)^ , 

where the loop functions are defined as 

fi{x,y) = ^j{i,x,y) + ^k{l,x,y), 

7 1 

f4{x,y) = -j{l,x,y) - -k{l,x,y), 

1 5 

fbix^y) = -j{l,x,y) + -k{l,x,y), (9) 
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and the j and k are defined in [18] . Tlie RG running of tfie Wiison coefficients down to nib 
scale can be found, for example, in [19]. 

We can calculate the Bg — Bg mixing matrix element, which is in the form 



SM, 



1 + R). 



(10) 



The mass difference of Bg — Bg system is then given by 



Arug = 2|Mf2| 



Am!^|l + i?|. 



(11) 



(12) 



In the SM contribution [20] to the mass matrix element 

127r2 

the non-perturbative parameters fs,, and Bb^ give main contribution to the theoretical 
uncertainty. Using the combined lattice result [21] from JLQCD [22] and HPQCD [23], 



f 



= (0.295 ± 0.036) GeV, 



(13) 



(HP+JL)QCD 



the SM predicts 

Amf^ = (22.5 ± 5.5) ps'^ (14) 

which is consistent with the values in (2) obtained from global fits. For the prediction in 

(14), we used riB = 0.551, mf^{mt) = 162.3 GeV and Vts = 0.04113 [24]. 

Now, inserting the CDF data in (1) and the SM prediction in (14) into (11), we obtain 

\l + R\ = 0.77i[J:J^^(exp) ± 0.19(th), (15) 

where the experimental and theoretical errors were explicitly written. The expression for 
R in our scenario is given by ^ 

RW = 6(/^6) +Cl(/^6) + T- 



4:Bi{nb) \mb{nb) + mg{nb) 



^The Bb, in (12) is related to -Bi(/Xb) as [20] 

Bb, ^ Bi(M6)[ai^HM6)]-'^/^' 



(16) 
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A Bi{nb) \mb{nb) + rusinb) J ^' 
where we defined (i = 1, • • • , 5) 

CM ^ cf^'\f.b)/cr{i^,), 

Uy^b) ^ cf^^{iib)/cf'{iib). (18) 

The relevant B-parameters are given in [25] by 

B.iiib) = 0.86(2) (It) , B^{^ib) = 1.17(2) , = 1-94(3) (If) . (19) 

Now we briefly discuss B Xsj constraint. The SUSY parameters we consider are 
also directly constrained by the measured branching ratio of inclusive radiative S-meson 
decay, B — > ^57. We take this constraint into account, although it is not expected to be 
so severe as in a scenario with LR or RL mixing. In the operator basis given in [26], the 
SUSY contributions to the Wilson coefficients of magnetic operators in our scenario are 



rSUSY _ 4 1 7ra^in2^ie^ t (r M 
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^susY _ 1 TTQis sin 26'Le"'^ 



{^\h{.Xb^g) - \jl{Xb^9)^ - {bL ^ Sl) 



At v^GFml 
where = V^^Vtb and 

1 — 6x + + 2x^ — 6x^ log X 



(20) 



h(x) = 



12(1 -X 



|4 



^, , 2 + Sx — 6x^ + + 6a; log a; 

Jlix) — : 71 . (21 

^ ' 12(1 -x)4 ^ ^ 

There are also chirality flipped Cj^^sg with L replaced by R. Therefore, we can see that in 
principle 9l(^r),Sl{K) and rris—m^ can be constrained. Compared to the LR{RL) mixing case 
where large SUSY contribution 0{mg/mb) is possible due to the chirality flipping inside the 
loop, our scenario allows only a small SUSY correction to the SM contributions. In addition, 
although LL mixing gives a hnear correction (^(C'f^fJ/C'l^sg) interference term, 

RR mixing generates only a quadratic correction C'dC'l^fJ/Cf^g^p) because it is added 
incoherently to the SM contribution. 



3 Numerical analysis 



In this Section, we perform numerical analysis and show the constraints imposed by Amf^^. 
We also consider the BR{B Xg'y) constraint. 




(a) (b) 

Figure 1: Contour plots for |1 + -R| in {nis^fiL) plane. Sky blue region represents 2cr allowed 
region (0.39 < |1 + i?| < 1.15), blue la allowed region (0.58 < |1 + i?| < 0.96), and white 
(grey) region is excluded at 95% CL by giving too small (large) Am^. The labeled thick 
lines represent the constant (bB}°\B Xsi) - BR^^{B X,-f)^ /BR^^{B X,7) 
contours. Only LL mixing is assumed to exist. The fixed parameters are rrig = 0.5 (TeV), 
m^^ = 0.5 (TeV), (a) 6l=0, (b) 6l = n/2. 

From (8) it is obvious that the larger the mass splitting between 's and b, the larger 
the SUSY contributions are. Therefore we expect that (15) constrains the mass splitting 
when the mixing angle O^pt) is large. This can be seen in Figure 1 where we show filled 
contour plots for |1 + i?| in (7715-^,6'^) plane: sky blue region represents 2cr allowed region 
(0.39 < \l + R\ < 1.15), blue la allowed region (0.58 < |1 + i?| < 0.96), and white (grey) 
region is excluded at 95% CL by giving too small (large) Arris. For these plots we assumed 
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that only LL mixing exists and fixed nig = 0.5 TeV, m^^ = 0.5 TeV. In Figure 1(a), we 
fixed Sl — 0- We can see that the SUSY interferes with the SM contribution constructively 
{i.e. the SUSY contribution has the same sign with the SM), and when the mixing angle is 
maximal, i.e. 9l — ±7r/4, m^^ ""^6^ cannot be greater than about 150 GeV. In Figure 1(b), 
we set 5l = 7r/2. The SUSY contribution can interfere destructively {i.e. in opposite sign) 
with the SM and much larger mass splitting is allowed. Therefore we can see that the 
allowed parameters are sensitive to the CPV phase. 

Also the constant (^BR^\B X,^) - BR^^{B X,-f)^ /BR^^{B X,-f) fines 
are shown. For fixed 9l, larger mass splitting mj^ — m^^ gives larger deviation for the 
branching ratio. This can be understood from (20). However, for very large mass splitting 
the SUSY contribution decouples and the deviation eventually saturates. We can see that 
BR^°^{B Xg^) deviates from the SM predictions at most about 5% in the region allowed 
by Am^. Since 

BR^''P{B X,-f)/BR^^{B X,-f) = 1.06 ± 0.13 (22) 

for > 1.6 GeV [27], it is clear that the BR{B Xgj) constraint is completely irrelevant 
in Figure 1. 

The plots for the scenario with RR mixing only are the same with Figure 1 because the 
expression for Bg — Bg is completely symmetric under L ^ R. As mentioned above, the 
contribution to BR{B Xgj) is much smaller than LL case. 

In Figure 2, contour plots for constant |1 + i?| in {9l,Sl) plane are shown. For Fig- 
ure 2(a)(2(b)), we fixed m^^ = 0.8(1.0) TeV. The other parameters used are the same with 
those in Figure 1. We can again see the strong dependence on the CPV phase 6l. It can also 
be seen that the parameter space with large mixing angle 9l can be made consistent with 
the experiments by cancellation with the SM contributions in the destructive interference 
region {i.e. 5l ~ 7''/2). 

Now we consider a scenario with both LL and RR mixing at the same time. Then the 
operators O4 and O5 are additionally generated as mentioned above. They dominate Oi 
or Oi for sizable mixing angles. As a consequence, the constraint on the SUSY parameter 



8 




(a) 



(b) 



Figure 2: Contour plots for |1 + R\ in {9l,Sl) plane, (a) mj^ = 0.8 (TeV), (b) m^^ = 1.0 
(TeV). The rest is the same with Figure 1. 



space is very stringent as can be seen in Figure 3. In Figure 3 we set rrig = 0.5 TeV, 
m^^ = m^^ = 0.5 TeV, m^^ = m^^ = 0.6 TeV, and (a) Sl = Sr = (b) 6l = 0,6r = tc/2. 
Even for small mass splitting most region of the parameter space is ruled out by giving too 
large Arris- We can see that BR{B — > Xg'j) is almost insensitive to the change of 6fi as 
mentioned before. 



4 The predictions of S^s and A 



^SL 



The CPV phase in the Bg — Bs mixing amplitude will be measured at the LHC in the near 
future through the time- dependent CP asymmetry 



(23) 
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-0.75-0.5-0.25 0.25 0.5 0.75 -0.75-0.5-0.25 0.25 0.5 0.75 

eL Ql 

(a) (b) 

Figure 3: Contour plots for |1 + B\ in {9l,0r) plane, m^^ = mj^ = 0.6 (TeV). (a) 6l = 
= (b) 6l = 0,6r = it/2. We assume both LL and RR mixing exist. The rest is the 
same with Figure 1. 

In the SM, S^^ is predicted to be very small, S^f = -sin 2/?, = 0.038 ± 0.003 = 
aj:g[(y^*Vtb) / (y*gVcb)]) [7]. If the NP has additional CPV phases, however, the prediction 

S^^ = - sin(2/3, + arg(l + R)) (24) 

can be significantly different from the SM prediction. 

In Figure 4, we show |1 + i?| constraint and the prediction of S^^ in (ms^,5i) plane. 
However, the B —>■ Xs'~f prediction is not shown from now on because it is irrelevant as 
mentioned above. For Figure 4(a), we assumed the scenario with LL mixing only and 
maximal mixing 9l = tt/A. We fixed m-g = 0.5 TeV, m^^ = 0.5 TeV. For Figure 4(b), we 
allowed both LL and RR mixing simultaneously, while fixing rrig = 2 TeV, m^^ = m^^ = 1 
TeV, rris^ = 1.1 TeV, 9r = tc/A, 5l = 7r/4, and 6ji = 7r/2. In both cases we can see that 
large S^^j, is allowed for large mass splitting between m^^ and mj^. At the moment, S^tf> 
can take any value in the range [—1,1] even after imposing the current Amf^^ constraint. 
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Figure 4: Contour plots for |1 + R\ in (ms^,^^,) plane. The 5*^^ predictions are also 
shown as thick contour lines. The thin red lines are constant y4|^[10^^] contours as- 
suming Re(r^2/^i2)^^ ~ —0.0040. (a) Only LL mixing is assumed to exist. We fixed 
m-g = m^^ = 0.5 TeV, 5^ = 7r/4. (b) Both LL and RR mixing are assumed to exist simulta- 
neously. We fixed m-g = 2 TeV, m^^ = m^^ = 1 TeV, m^^ = 1.1 TeV, 9r = 7r/4, 6l = vr/4, 
and 6r = tt/2. The rest is the same with Figure 1. 

Finally we consider the semileptonic CP asymmetry [28, 16, 7] 
= r(g. ^ £+X) - TjB, i-X) ( VI, \ 



It is approximated to be [7] 

where Re(r^i2/^f2)^^ = -0.0040 ± 0.0016 [29]. The SM prediction is A|l(SM) = (2.1 ± 
0.4) X 10-5 [29, 30]. 

In Figure 4, the thin red lines are constant ^5l[10-^] contours taking Re(ri2/Afi2)^^ = 
—0.0040. We can readily see that the strong correlation between S^^ and Agj^. This can 

II 




Figure 5: The correlation between and S^^. The red hne is 1-cr upper bound, 
be seen from the relation 



4* 



Re 
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Ml, 



(27) 



\l + R\ 

For small R the two observables are linearly correlated as can be seen in Figure 4. 

In Figure 5, we show the correlation between and 5*^^. We scanned 0.5 < m-g < 4.0 
TeV, 0.5 < m~^^,m~s^ < 2.0 TeV, -7r/4 < Ol < 7r/4 and < 6l < 2n, while fixing 
rrig = m^^ = 0.5 TeV. The Am^ constraint is imposed with 0.39 < |1 + -R| < 1.15. We have 
checked that in the scenario with only LL (RR) mixing, we get the similar correlations. 
The red line is experimental 1-a upper bound from = —0.013 ± 0.015 [7]. Now several 
comments are in order: i) The values for S^^ and can be significantly different from 
the SM predictions, ii) The two observables are strongly correlated. These two facts were 
already noted in [7] . It has been checked that in the {KeR, ImR) plane the above scanned 
points can completely fill the region allowed by Am^. This explains why the correlation in 
Figure 5 is basically the same with model- independent prediction in [7]. iii) Although it 



looks like that large negative S^^i, value is disfavored, due to large error in Re{Tl2/M 
we cannot definitely rule out the region at the moment. 
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5 Conclusions 



We considered the MSSM scenario with large LL and/or RR mixing in the down-type 
mass squared matrix. This scenario is strongly constrained by the recent mesurements of 
Bg — Bg mass difference, Amg, in contrast with the MSSM scenario where the flavor mixing 
is controlled only by the CKM matrix [16, 8]. The constraint is most stringent when both 
LL and RR mixing exist simultaneously. It is also shown that the allowed region is quite 
sensitive to the CP violating phase. 

We also considered the time-dependent CP asymmetry, S^^, and the semileptonic CP 
asymmetry, A^^^. It was shown that the S^^ and A^j^ can take values significantly different 
from the SM predictions. There is also strong correlation between S^^ and AgL. 
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